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Abstract 


The report discusses momentum management for the CDG Planar Space Plat- 
form. The external torques on the Space Station are assumed to be gravity 
gradient and aerodynamic with both having bias and cyclic terms. The inte- 
grals of the cyclic torques are the cyclic momenti which will be stored in 
the momentum storage actuator. Various techniques to counteract the bias 
torques and center the cyclic momentum were investigated Including gravity 
gradient desaturation by adjusting vehicle attitude, aerodynamic desatura- 
tion using solar panels and radiators and the deployment of flat plates at 
the end of long booms generating aerodynamic torques. 


Nomenclature 


A p transformation from p-th appendage principal axis to core body frame 

cm center of mass 

C d drag coefficient (assumed = 2) 

CMG control moment gyro 

cp center of pressure 

GG gravity gradient 

H instantaneous stored momentum 

jT estimated average angular momentum 

I identity matrix 

I0P in-orbit plane 

Ip inertia of p-th appendage 

I pcm inertia of p-th appendage in cm coordinates 

Iy total vehicle inertia tensor 

K d p atmosphere diurnal bulge 

LV local vertical 

nip mass of p-th appendage 

fl unit normal to p-th appendage 

N-m newton meters 

N-m-s newton meters seconds 

01 orbit inertial 

POP perpendicular to orbit plane 

£ unit radius vector from center of earth to center of mass 

S surface area of p-th flat plate 

T B y desired y-axis torque for bias counteraction 

Tq desired control torque 

average value of I0P torque in desired direction 
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Nomenclature (Continued) 

Tqiop d ® s ^ re d control torque magnitude In orbit plane 

T gg gravity gradient torque 

Tjy desired y-axis centering torque 

T n0|T) nominal torque from one appendage if p « p nom and air flow is per- 
pendicular to surface 

T 0 orbital period 

Tp a aerodynamic torque due to p-th appendage 

T average value of IOP torque perpendicular to desired direction 
U vehicle angle from orbital noon 

peak atmospheric density phase lag 
U on half angle that IOP aerodynamic desaturation is "ON" 

U T angle of desired torque IOP with orbital noon 
V velocity magnitude 

9 unit vector along relative wind 

Xp moment arm of p-th appendage 

z z-transform variable 

Al^ inertia difference (aI^I^-Ij, i, j, k in even permutation) 

a«j» ^ A-rotation of vehicle about POP axis to remove bias 

A<t >2 A-rotation of vehicle about POP axis to center momentum 

a angle at panel "ON" makes with air flow for POP desaturation 

6 angle of sunline to orbit plane 

y angle of "ON" POP desaturation panel with air flow 

<t> angle of vehicle z-axis to local vertical about POP axis 

<fip angle of vehicle principal z-axis to local vertical about POP axis 

ui 0 vehicle angular rate in orbit 

it 3.1416 
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Nomenclature (Continued) 

0 instantaneous atmospheric density 

0 avg averagp atmospheric density 

o est estimated average atmospheric density 

p nom nominal atmospheric density 

o vehicle x-axis rotation for GG desaturation 

o' maximum o adjusted for finite maneuvering rate 

®max max rate for maneuverin 9 vehicle about its x-axis 
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1.0 Introduction 


Figure 1 represents the CDG planar space station configuration whose basic 
attitude has the Z-axis local vertical (LV), the Y-axis perpendicular to 
the orbit plane (POP) and the X-axis along the velocity vector. The solar 
arrays are maintained nominally perpendicular to the sun? me requiring two 
degrees of freedom, one for the vehicle rotation and the second for the 
sunline with respect to the orbital plane ( 0 -angle), The radiators are 
maintained POP and only edgewise illuminated by the sunline. However, un- 
like the solar arrays which have unlimited rotation capability, due to the 
necessity of fluid transfer between the space station and the radiators and 
the limited freedom of the plumbing interface, the radiators must be un- 
wound during the dark side of the orbit. 

If the balanced vehicle defined In Figure 1 were ideal, the principal axes 
of the core body would be perpendicular to the orbit plane and along local 
vertical and no gravity gradient torques due to core body would exist. The 
center of mass would be at the geometric center and no aerodynamic torques 
would exist. The only remaining torques would be a bias torque about the 
X-axis due to the non-zero e-angle of the solar arrays and a cyclic torque 
about the POP axis due to the solar arrays. The radiators only contribute 
torque during the unwind cycle and were determined not to be significant 
for sizing momentum storage and desaturation requirements. For the actual 
inertias and center of mass used for the study, significant core body grav- 
ity gradient biases and aerodynamic torques do exist although the balanced 
design insures their being of small magnitude for the size of the vehicle 
investigated. The momentum storage system assumed for this study was four 
Skylab type double gimbal control moment gyros (CMGs) having an angular 
momentum of 3000 N-m-s each. 

The main effort of this study was to develop momentum management techniques 
to counteract the bias torques and also to center the stored cyclic momen- 
tum along three inertial axes. Prior to developing momentum management 
techniques, it was necessary to define disturbance torques for the basic 
vehicle. Two types of disturbance torques were admitted for this study; 
gravity gradient and aerodynamic. Gravity gradient torques were defined by 
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FIGURE 1 CDG PLANAR SPACE STATION 
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the standard equations while aerodynamic torques were defined by a combina- 
tion of flat-plate areas, a simple orbital curve fit to the Jacchia density 
model and Inelastle collision of particles with the flat plates. The defi- 
nltion of torques Is complicated by the relative motion of the appendages 
with respect to the core body leading to the selection of seven independent 
bodies to represent the total Inertia. These bodies are the core, four 
solar arrays and two radiator sets. Nominally, the four solar arrays and 
two radiator sets each move as units but separate motions were allowed In 
order to investigate their use as desaturation techniques. 

2.0 Description of External Torques 

The inertia of the entire vehicle Is defined as that of the core body plus 
the inertias of the appendages reflected to the center of mass of the sys- 
tem. Due to their symmetric nature, the center of mass of the solar arrays 
and the radiators are at the center of symmetry of the core body. However 
the ?r of mass (cm) of the core body is slightly offset from the center 
of ymmetry. The inertia of an appendage relative to the cm of the core 
body consists of the sum of the inertia due to its own inertia matrix about 
its cm and the Inertia moment of a point mass at the appendage cm acting 
about the core body cm. The appendage inertia must be transformed via a 
similarity transform from its principal axes to the axes fixed in the core 
body. The transformation varies with appendage attitude. The variable 
inertia due to an appendage is defined by equation 1 while the gravity gra- 
dient torque is defined by equation 2. The gravity gradient torque on the 
total vehicle can be defined two ways of which the first is the gravity 
gradient torque due to the total inertia of the core body and appendages. 
The second is to define the inertia of each appendage, determine the gravi- 
ty gradient torque and add the individual torques to generate the total 
gravity gradient torque. The latter approach allows one to see the contri- 
bution of the individual appendages for insight Into possible desaturation 
techniques. The aerodynamic torque due to the core body and each appendage 
is simply the force on each body and a vector cross product of that force 
with the vector between the cm and the center of pressure (cp). This 
torque is defined by equation 3. 
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The flat plate area of the core is defined as the projection of a cylinder 
onto a plane. Figure 2 shows the aero torques due to the core body. The 
torque perpendicular to the orbit plane (POP) is in the same direction and 
only varies in magnitude due to the diurnal bulge. The torques in the or- 
bit plane (IOP) are essentially cyclic but have bias components due to the 
diurnal bulge. The gravity gradient (GG) POP torque shown in Figure 3 is a 
bias while the IOP torques are cyclic in Orbit Inertial (01) space. All 
IOP torques and momenta will be displayed in 01 rather than body axis as 
that Is the frame where momentum buildup truly occurs and in which it must 
be desaturated. The density variation in an orbit is given by equation 4, 
a simple curve fit to the Jacchia model of the nominal atmosphere. 


p = p avg (1 + f 'DB cos ( u " U DB^ 


( 4 ) 


Figure 4 shows a typical density variation over an orbit with the peak den- 
sity occurring 30 degrees after noon and the minimum density occurring 30 
degrees after midnight. Figure 5 shows the aero torque due to a single so- 
lar array with significant torques on all axes which is the same condition 
observed for the GG torques of Figure 6. Figure 7 shows that the aerody- 
namic torque due to four solar arrays is much less as the majority of the 
moment arms cancel. If the cm of the core body were at the center of sym- 
metry, no solar array aero torques would exist. The cyclic gravity gradi- 
ent torques due to the individual panels cancel IOP but add POP. Figures 
8, 9, 10 and 11 show that the total torques on the vehicle are dominated by 
the solar arrays. However the dominant bias component POP is the GG torque 
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FIGURE 7 AERO TORQUE, FOUR SOLAR ARRAYS 
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. FIGURE 8 GG TORQUE, FOUR SOLAR ARRAYS 
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FIGURE 11 TOTAL TORQUE, TOTAL VEHICLE 
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due to the core body. The dominant IOP bias contribution is the solar ar- 
rays aero torque at a e-angle of zero while solar array 6G torque dominates 
at e-angles of magnitude close to 45 degrees. Figure 12 shows the stored 
momentum profile if no momentum management were active and the Initial mo- 
menta were -1000 N-m-s In Orbit Inertial (01) X and l respectively and 
-3600 N-m-s along the Y-axis which is perpendicular to the orbit plane 
(POP). The 01 X-axis Is along the velocity vector at noon and the 01 /-ax- 
is points up at noon. The start of all runs is orbital noon with the core 
body approximately trimmed to remove the POP bias torque. 



FIGURE 12 MOMENTUM WITH NO MOMENTUM MANAGEMENT 
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3.0 Techniques Considered 


A known technique for bias counteraction is the use of magnetic torquing 
using four Space Telescope magnetic torquer bars. These wc Id have the 
capability of generating approximately .13 N-m torque per axis when aver- 
aged over an orbit and are Independent of vehicle configuration, att/tude 
and atmospheric density. As their capabilities and control laws are known, 
including the disadvantages of weight, power and magnetic contamination, 
their use was not specifically Investigated in this study but they remained 
the benchmark for innovative techniques. Several innovative techniques 
were investigated with various degrees of success. The emphasis of the 
study was to Investigate in depth those techniques that appear promising. 
Among techniques that were investigated but eliminated were using the in- 
wind profile of the radiator which did not generate sufficient torque. 
Modifying the pointing of various solar arrays using both the a-angle gim- 
bal and the local vertical gimbal did not produce adequate desaturation 
without very large rotations during the dart side of the orbit which was 
considered undesirable from 'inting performance. Off-nominal pointing of 
some of the solar arrays dr ,ig daylight would produce undesirable power 
profiles. 

The promising techniques developed to be described are use of the core body 
and the deployment of special aerodynamic appendages. All the techniques 
developed would be autonomous in the flight computer and simple to imple- 
ment requiring only knowledge of inertias and approximate average atmos- 
pheric densities. Since the core body is maintained nominally local verti- 
cal (LV), its use for POP desaturation by means of gravity gradient torques 
is well known. The control laws for this will be described. Phased mo- 
tions of the core body about the vehicle X-axis can be used to desaturate 
momentum IOP and will also be described. Initially, the control concept 
was to make these desaturation maneuvers every orbit. Results showed that 
after the initial transient decayed, the maneuvers were very small although 
payload pointing interference was still possible. Therefore an alternate 
control concept was to allow absolute momentum to approach a value inside 
its capability sphere. Only at this point would desaturation maneuvers of 
the core body be required although they would be much larger than those 
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performed each orbit. It appears that single orbit maneu /ors performed 
every tenth orbit might be adequate after initial equilibrium is attained. 


Using the appendages of Figure 13, no core body motion is required. Closed 
loop autonomous control laws for the angle of attack of <hese appendages 
have been developed. The appendages along the vehicle Y axis desaturate 
the I0P momentum while the appendages along the Z-axis control the POP mo- 
mentum. The hardware required to deploy and control these appendages could 
be adapted from solar array and antenna drives and does not represent new 
technology. The attainment of a given torque at a certain atmospheric den- 
sity depends on the product of area times moment arm which could be traded 
off depending on considerations of weight, complexity and possible inter- 
ference with payload pointing. The results of using these appendages com- 
pares favorably with magnetic desaturation in terms of weight, power, and 
complexity of both hardware and control laws. When atmospheric density is 
high, the baseline sized appendages produce greater desaturation capability 
than four Space Telescope torquer bars. During period of minimum density, 
the torquer bars have more capability but the requirements are less. 
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IOP DESMURATION 
APPENDAGE 1 
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POP DESATURATION 
APPENDAGE 2 


CDG PLANAR SPACE STATION WITH AERODYNAMIC DESATURATION APPENDAGES 
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4.0 Gravity Gradient Desaturation Using the Core Body 


Utilizing the core body, which maintains an attitude relatively close to 
local vertical (IV) to desaturate the momentum perpendicular to the orbit 
plane (POP), is a technique previously investigated (ref 1). Assuming one 
principal axis is approximately POP and a second is approximately LV, the 
gravity gradient torque along the POP axis is given by equation 5. The 
total desired attitude change is the sum of that required to remove the 
bias and that required to center the momentum along the POP axis. The at- 
titude change to counteract the bias torque is given by equation 6 which 
also utilizes the inverse of equation 5. The angle change required to cen- 
ter the momentum is given by equation 7. The new attitude about the POP 
axis is given by equation 8 and includes an effect for increased momentum 
due to the bias acting over an additional half orbit. This angle change is 
accumulated in the flight computer and while the core body attitude about 
the POP axis is controlled within a small range by the requirements of 
momentum management, the flight computer will, at all times, have knowledge 
of the actual attitude. Figure 14 shows a typical attitude about the POP 
axis while Figure 15 shows the momentum. The large initial changes about 
the POP axis occur due to the requirements of the centering logic. 


Tqgy = 


(5) 
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FIGURE 15 POP MOMENTUM WITH CONTINUOUS CORE BODY CONTROL 
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If the vehicle Z-principal axis Is moved out of the orbit plane by a rota- 
tion about the vehicle X-axis* a body fixed gravity gradient torque about 
the X-axis will be produced. However this torque is cyclic about two axes 
fixed in the orbit plane. In order to generate a gravity gradient bias 
torque in the orbit plane, it is necessary to generate a phased rotation 
about the vehicle X-axis. If the orbit inertial (01) axes have average 
momenta over the previous orbit given respectively by H x and W z the total 
bias torque to be generated about the 01 X and Z axes are given by equa- 
tions 9 and 10. The magnitude of the total IOP torque is given by equation 
11 and the direction of the desired torque is given by equation 12. 



( 10 ) 


T, 


DIOP 


_ / 


'Dx 


■+ T 


Dz 


( 11 ) 


Uj = tan 


-1 



( 12 ) 


After the desired magnitude and direction of the IOP torque to be generated 
is calculated, the second part of the calculation is the definition of the 
vehicle attitude profile with respect to the orbital plane. Instantaneous- 
ly, the gravity gradient torque in body space is given by equation 5 with 
All replacing aI 2 . If that attitude is maintained over one-half orbit cen- 
tered about a given direction, the average torque in that direction is giv- 
en by equation 13 while the average torque in the perpendicular direction 
(equation 14) is zero. Reversing the rotation about the X axis during the 
other half of the orbit produces a rectified torque in the desired direc- 
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tion and still zero in the orthogonal direction. The angle required to 
desaturate the I0P bias is given by equation 15 which shows an efficiency 
less of 2/tt compared with the angle required for POP gravity gradient 
desaturation. Examining Figure 16, which shows the IOP momentum with core 
body desaturation active, indicates that the average values of 01 X and Z 
momenta impose only small momentum storage requirements on the CMGs and 
that the bias torques are removed. Figure 17 shows a typical core body 
rotation profile about its X-axis for IOP desaturation which shows that 
after the transients have decayed, the required angles are quite small and 
tiiet even maneuvering at .01 deg/sec, the times to attain the new attitudes 
are negligible. 
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It should be noted that since torque and momentum are not in the same 
direction, it is difficult to remove both bias and average values of the 
momentum in two 01 axis. However the average momenta are at values close 
to zero and in maintaining those values with the control law, bias torques 
are also counteracted. Examining the POP rotations of Figure 14 and the 
IOP profile of Figure 17 shows that after the momentum along POP is cen- 
tered and the bias is removed, both the POP corrections required and the 
IOP profiles represent small angles which, nevertheless, could disturb pay- 
load pointing. Due to the balanced structure of the vehicle, the momentum 
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accumulation, on a per orbit basis, is very low once POP bias h-s been re- 
moved by the appropriate rotation of the core body. This led to the de- 
velopment of an alternate method of gravity gradient desaturation using the 
core body. Instead of performing an adjustment every orbit, the vehicle 
was held inertial ly fixed until the absolute momentum attained a software 
limit. The maximum momentum storage capability of four Skylab type Control 
Moment Gyro (CMG) is 12000 newton-meter-seconds (N-m-s) leading to the se- 
lection of a software limit set at 9000 N-m-s. The logic was set so that 
if the stored absolute momentum attained this value, corrective action for 
both PGP and IOP momentum would be performed. 

The amount of momentum left for maneuvering is quite low so that a maximum 
core body maneuvering rate of .01 deg/sec about any vehicle axis v/as as- 
sumed. Since the angles required for desaturation will be much larger when 
momentum is allowed to go to the inside of the 9000 N-m-s momentum sphere, 
the assumption of square wave attitude angle profiles of Figure 16 for IOP 
desaturation can no longer be made. Therefore with the assumed maximum 
core angle with infinite rate capability calculated, the actual angle is 
adjusted for finite rate according to equation 16. This maintains the same 
integral of area x time as if the new attitude were attained instantaneous- 
ly leading to the trapezodial profile of Figure 18. 

o' - V xl W (l.^. 80 J (16) 

°max^o/2^ 

These considerations are less critical for the new angle about the POP ax- 
is. Desaturation about the POP axis is not dependent on angle in orbit for 
the adjustment of the momentum states. Therefore finite maneuvering times 
will compensate at the beginning and end of the desaturation profile. Due 
to the significant time required for the completion of the single orbit IOP 
desaturation, it was decided to perform a maneuver in only two steps as 
shown in Figure 18 rather than the three steps of Figure 17. It should be 
noted that Figure 18 shows the desaturation profile IGP from time of start, 
which could be any time in orbit, and assuming the calculations are per- 
formed at a fixed position in orbit, the end of the desaturation could take 
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FIGURE 18 

TILT PROFILE FOR SINGLE ORBIT DESATURATION, RATE LIMIT = .01 DEG/SEC 









two orbits from the time that the calculations are made. The other change 
in IOP desaturation when no longer performed at orbital frequency is, rath- 
er than drive the IOP momentum to zero when desaturation is required, the 
momentum will be driven to the opposite point on the IOP momentum circle 
given by equations 17 and 18. Since the bias should drive the momentum in 
the same direction as previously, the maximum time allowed before the next 
desaturation is required should be extended. Figure 19 shows the momentum 
IOP about two inertial axis with single orbit desaturation. 



(17) 
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Dz 



(18) 


The single orbit POP desaturation technique utilizes only minor modifica- 
tions to the POP desaturation technique at orbital frequency. The angle to 
remove the bias is given by equation 6 and that remains as the new vehicle 
attitude until the next desaturation is required. The angle to center the 
momentum is given by equacion 7 and that angle is removed at the conclusion 
of the single orbit desaturation. The new angle at the beginning of a de- 
saturation is given by equation 19 and that at the end of one desaturation 
period is given by equation 20. The period of the single orbit POP desatu- 
ration is concurrent with that of IOP desaturation which generally is not 
the position (orbit noon) where the calculations are performed. 


4> — 1 • 5a<(>^ - A4>2 
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The POP attitude maneuver at the beginning of the single orbit desaturation 
is the sum of the bias removal and momentum centering angles while that at 
the conclusion is the negative of centering angle. Using the momentum de- 
saturation technique of allowing absolute momentum to reach the maximum al 
lowable value of 9000 N-m-s, the POP momentum with single orbit desatura- 
tion is shown in Figure 20. A subsequent desaturation was not required for 
30 orbits which is probably too optimistic for a real mission as simulation 
did not include varying e-angle, atmosphere density and mass properties 
between successive orbits. The real trimming of torque and momentum will 
probably not be as accurate as shown in the simulation. However even the 
attainment of ten orbits of pure Inertial Hold before momentum management 
trim maneuvers might represent adequate viewing time for many payloads. 

Another desaturalion technique, magnetic desaturation using the control 
laws developed during the Space Platform study was considered but not actu- 
ally investigated in detail for Space Station. These control laws would be 
directly applicable to Space Station and the magnitudes of the I0P bias 
torques (.05 N-m) on 01 X and Z are well within the capability of four 
Space Telescope torquer bars. However the .5 N-m bias torque about the POP 
axis is beyond the capability of 4 bars unless the vehicle is first trimmed 
by a closed loop POP axis rotation. Magnetic desaturation also has unde 
sirable properties such as the weight and power required and the potential 
for magnetic contamination of payloads. 

Two alternate desaturation techniques would be mass shifting for gravity 
gradient desaturation and appendage control for aerodynamic desaturation 
with the latter considered more promising. The appendages considered for 
this study are two flat plates along ± vehicle Y-axis for I0P desaturation 
and two flat plates along ± vehicle Z-axis for POP momentum trimming. It 
was decided that the magnitude of the POP bias was too large for complete 
aerodynamic desaturation and the flat plates would act as trim tabs after 
several orbits of gravity gradient maneuver desaturation about the POP ax- 
is. This report will discuss sizing, control laws and show closed loop re- 
sults of aerodynamic desaturation using flat plate appendages. In addition, 
results will be shown for those cases where the atmosphere is different 
from the nominal density used to size the appendages and control law gains. 
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5.0 Aerodynamic Oesaturation Using Flat Plate Appendages 


Simulation results have indicated that for the baseline vehicle inertias, 
baseline offsets of the center of mass from the center of symmetry and the 
baseline atmosphere model (o aV g s 1*73 * 10”* 2 kg/m 3 , Kdb = .h, Udb 30 
degrees), the 01 X and Z bias torques are about .0135 N-m each and the mag 
nitude of the total IOP bias torque which is the vector sum of the compo- 
nents is about .02 N-m at a zero degree B-angle. However at a non-zero 
angle, although the aerodynamic forces and torques are less, the gravity 
gradient bias IOP is greater. The worst case IOP bias torque occurred at a 
6-angle of -45 degrees about the X-axis. This can be realized intuitively 
that as the solar arrays are rotated from POP for the 0-angle, a gravity 
gradient bias torque along the 01 X-axis is generated having its maximum 
value at 45 degree magnitude s-angle. For the inertias and cm used, the 
magnitudes of the aerodynamic and gravity gradient biases add about the 01 
X-axis at the negative s-angle producing the worst case torque of .045 N-m 
consisting of 0.035 N-m gravity gradient and 0.01 N-m aerodynamic torque. 
Using the efficiency factor of 2/n (to be shown later) which is identical 
to that for IOP desaturation using the core body, the required torque from 
the panels would be about .075 N-m. To allow margin, the panels for IOP 
desaturation were sized to generate T nom of 0.1 N-m. Assuming an orbital 
velocity squared of 5.8 * 10 7 meters 2 /sec 2 and a drag coefficient (Cj) of 
two, the product of area times moment arm must equal 1000 meters cubed to 
generate T nom . Typically, this would represent an area of 40 square meters 
with a 25 meter moment arm or 20 square meters with a 50 meter moment arm. 

The sizing requirements for the POP axis trim tabs is less clear-cut as re- 
quirements depend on the degree of trim attained by gravity gradient desat- 
uration using the core body. In order to simplify hardware considerations, 
it was decided to use the same requirements as for IOP desaturation, namely 
r nom torque capability at an atmospheric density of o nom . In the POP axis 


35 


aerodynamic torque, the unfavorable factor of 2/n for the theoretical 
torque achieved in the desired direction does not ^xist so the 'orque at- 
tained could be l nom - Referring to Figure 13, panel 2 generates a positive 
torque about the POP axis. The general requirements for deployment and ro- 
tation of the panels are consistent with solar array and antenna mecha- 
nisms. Candidates for actuators Include stepping motors and torquers. 
Since the vehicle is rotating at orbital rate, in order to maintain a fixed 
orientation to the air flow, the panels must only be rotated when reorien- 
tation is required. This could favor stepper motors as actuators since the 
devices would be off when not actually rotating. An estimate of the weight 
is 15kg for each actuator/boom/panel or 60 kg for the system. 

The control law for IOP desaturation was based on assuming a constant at- 
mospheric density and sampling the momentum of the two IOP axis at orbital 
frequency with the sampling performed at orbital noon. The choice of noon 
for momentum sampling is arbitrary for this simulation only. Perhaps mid- 
night which was used in Skylab for Ephemeris updating would be operational- 
ly preferable but results should be independent of sampler placement in or 
bit. The requirement of knowledge of average density makes the aerodynamic 
desaturation somewhat less desirable than gravity gradient desaturation. 
Currently, the inertia matrix of a spacecraft is assumed to be well known 
while atmospheric data at the Space Station altitude is known best from its 
long term averages. Therefore the study involving aerodynamic desaturation 
with flat appendages will consider the case where the control law estimate 
pi o aV q is different from the "real" atmospheric density. It will be shown 
that momentum desaturation is still viable even with a certain amount of 
mismatch. The density variations within the orbit, due to the diurnal 
bulge are less critical than the average value of the density for aerody- 
namic desaturation and because of the uncertainty of the exact density dis- 
tribution during an orbit, are not included in the control law. 
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The momentum sampling scheme for the generation of the required torques 
follows along the same general lines as those used for the core body desat- 
uration with the IOP bias torques to be counteracted given by equations 9 
and 10, the magnitude of the total bias t^ue to be generated is given by 
equation 11 while the direction of the d^.red torque is given by equation 
12 . 


Figure 21 shows the two appendages for IOP desaturation with 2 "on" and 1 
feathered producing a negative torque about the vehicle Z-axis. When ap- 
pendage 1, along vehicle +Y-axis is rotated perpendicular to the air flow, 
a positive torque about vehicle Z is generated. The appendages which de- 
velop torque along the vehicle Z-axis IOP are along the two ends of the ve- 
hicle Y-axis with their .'.is POP. The control law developed for IOP desat- 
uration admits on-y three states for the two appendages: both feathered, 
appendage 1 perpendicular and 2 feathered and 1 feathered and 2 perpendicu- 
lar th the velocity vector. Therefore only three torque states about Z-ve- 
hicle axis are possible. These are zero torque if both or none of the ap- 
pendages were feathered, and maximum magnitude if only one appendage is 
feathered. The zero torque condition will be assumed by both appendages 
feathered. If the same appendage were perpendicular to the air flow for 
the entire orbit, the torque would be a bias in vehicle space but cyclic in 
IOP axes. Therefore it will be necessary to develop a control law involv- 
ing feathering each panel part of the orbit and possibly both panels part 
of the orbit to generate a torque with desired magnitude and direction. 

It is assumed in this study that the separation of the cm from the cp is 
negligible compared to the moment arm of the individual appendage and given 
identical densities, the torques due to each appendage would have the same 
magnitude. It was also assumed for this study that the appendages do not 
affect the gravity gradient torques. A typical case that is also easy to 
visualize is a e-angle of zero where the bias torques are -.014 N-m about 
both 01 X and Z axes. This also represents a torque of magnitude .02 N-m 
about an axis 315 degrees from noon. In order to counteract this external 
bias, the appendages should generate a positive torque along an axis 45 
degrees from noon without any component along a perpendicular axis. 
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This is done by using panel 1 symmetrically about 315 degrees and panel 2 
about 45 degrees. All subsequent calculations are based only on density 
equal to p nom with closed loop calculations taking care of actual atmos- 
pheric density. The theoretical torques generated are symmetrical in 4 
quadrants about the torque direction so the equation for turning panels on 
is given only for two quadrants symmetric about the desired torque direc- 
tion. 
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FIGURE 21 IOP DESATURATION APPENDAGES 
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The average torque In the desired direction over one half orbit is given by 
equation 21. If a certain magnitude is desired, the inverse of equation 21 
is solved for the angle U Qn which represents the "on time" of the appendage 
torque I0P. This angle is given by equation 22. The maximum torque that 
can be generated in the desired direction if the angle is ninety degrees is 
T nom x 2 /h. The torque in the direction perpendicular to the desired 
torque is given be equation 23 and can be seen to be zero if p is actually 
constant. 
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i on 
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1 * p nom U T - 1 

T on 

The POP axis control contains an integr 1 path so the desired setting is 
the sum of a bias and integral path. The bias contribution to the desired 
POP torque is given by equation 24 while the desired centering torque is 
given by equation 25. The total torque desired is the sum of the bias and 
centering torques. The POP axis trim tabs control law utilizes a different 
approach from the IOP control law. Either panel 1 or panel 2 is "on" the 
entire orbit with the same angle based on the POP momentum sample. Unlike 
the IOP control law, the panel that is "on" can be set at any angle between 
zero and ninety degrees. The magnitude of the torque is determined by equa- 
tion 26 while the angle that the normal to the panel makes with the airflow 
to generate the desired torque is given by equation 27. 
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The sign of the desired torque determines which panel is "on". Figure 22 
shows the IOP appendage torque in vehicle space for a typical crbit, Figure 

23 shows the torque in the desired and perpendicular direction and Figure 

24 shows the torque in 01 axes. The torques shown are not perfect as they 
represent the torques generated by the control law using the first momentum 
samples. The closed loop nature of the calculations will assure that after 
every several orbits, t!-e profile converges to produce the torque in the 
desired direction despite the non uniformity of the atmosphere. This non 
uniformity can be seen in Figure 22 where a constant p would produce a 
square wave profile. 

Figures 25, 26 and 27 show IOP momentum buildup for nominal, 20% nominal 
and 5 times nominal average density over 9 orbits without bias counterac- 
tion. Figures 28, 29 and 30 show IOP momentum when the average density, 
nominal, 20% nominal and 5 times nominal, used in the control law exactly 
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FIGURE 23 DESIRED AND PERPENDICULAR IOP APPENDAGE TORQUE 
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FIGURE 25 

IOP MOMENTUM, NO APPENDAGE CONTROL, RHOAVG = 1.7E-12, BETA = -45 
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FIGURE 26 B 

IOP MOMENTUM, NO APPENDAGE CONTROL, RHOAVG = .34E-12, BETA = -45 { 
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FIGURE 28 

IOP MOMENTUM, APPENDAGE CONTROL, TRIMMED CORE, 
RHOAVG = 1.7E-12, RHOCNTRL = 1.7E-12, BETA = -45 


47 








48 







§ 

h 

FIGURE 30 

I0P MOMENTUM, APPENDAGE CONTROL, TRIMMED CORE, 

RHOAVG = 8.5E-12, RHOCNTRL = 8.5E-12, BETA = -45 
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matches the average atmospheric density. At the o=anglp of =45 degrees, 
f igu«"e ? 9 shows that regardless of the control law, the appendages can not 
counteract the I0P bias about the 01 X-axis bias torque when the atmospher- 
ic density is only 20% nominal. The only solutions are to increase the pan- 
el size or raise the solar arrays for 0 degree e-angle as shown in Figure 
31. This would reduce the available power if the actual e-angle were -45 
degrees. Figure 32 shows IOP momentum when the o es{ . assumed by the control 
law is different from the true o avq . The true density is 5 times that used 
by thn control law and makes the IOP bias desaturation marginally stable. 

T he performance of the panels for POP desaturation is heavily dependent on 
the approximate balance of the vehicle gravity gradient and aerodynamic bi- 
as torques. Figure 33 shows the approximate POP axis trim as a function of 
o aV g where vehicle trim is the angle between the vehicle Z-axis and the lo- 
cal vertical. Figure 34 shows the POP appendage torque when one is perpen- 
dicular to the air flow and with the nominal atmospheric model. The appen- 
dages are capable of generating .1 N-m torque at the nominal p aV g which 
represents trim capability over one degree. At 5 times the nominal densi- 
ty, the trim capability is 5 degrees while at 20% the trim capability is 
only .2 degrees of vehicle attitude. Therefore reaching an approximately 
balanced attitude about the POP axis using gravity gradient desaturation is 
most critical during periods low atmospheric density. Figures 35, 36 and 
37 show the POP momentum for an approximately trimmed vehicle at nominal, 
20% nominal and 5 times nominal density without additional appendage desat- 
uration. Approximately trimmed means the the POP bias is not perfectly bal- 
anced but within the capability of the aerodynamic desaturation panels when 
the density in the control law is exactly matched to the average density. 
Figures 38, 39 and 40 show typical POP momentum where the vehicle is 
trimmed and the true dens'ty and that assumed in the control law are 
matched. Figure 41 shows the POP momentum when the true p aV g is 5 times 
that assumed by the control law. Only a slight oscillation can be seen so 
that POP desaturation is quite insensitive to a p aV g greater than that as- 
sumed by the control law. However if the p aV g is less than that assumed by 
the control law, the degree of core body trim might be inadequate and the 
POP momentum would diverge. The divergence is not caused by lack of ro- 
bustness in the controller but by insufficient torque capability. 
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FIGURE 31 

{ IOP MOMENTUM, APPENDAGE CONTROL, TRIMMED CORE, 

RHOAVG = .34-12, RHOCNIRL = .34-12, BETA = 0 
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FIGURE 32 

IOP MOMENTUM, APPENDAGE CONTROL, 
RHOAVG = 8.5E-12, RHOCNTRL = 1.7E-12, BETA 












FIGURE 33 

POP AXIS TORQUE EQUILIBRIUM ATTITUDE AS A FUNCTION OF AVERAGE DENSITY 
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FIGURE 34 POP APPENDAGE TORQUE 
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FIGURE 35 

POP MOMENTUM, NO APPENDAGE CONTROL, TRIMMED CORE, RHOAVG = 1.7E-12 
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OP MOMENTUM, flPPENDOGE CONTROL, TRIMMED 
RHOfi VG= . 34E— 1 2 . RHOCNTRL= . 34E- 12 , BETft=- 



FIGURE 39 

POP MOMENTUM, APPENDAGE CONTROL, TRIMMED CORE, RHOAVG = .34E-1Z 
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FIGURE 40 

I0P MOMENTUM, APPENDAGE CONTROL, TRIMMED CORE, 
RHOAVG = 8.5E-12, RHOCNTRL =* 8.5E-12, BETA = -45 
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FIGURE 41 

POP MOMENTUM, APPENDAGE CONTROL, TRIMMED CORE, 
RHOAVG = 8.5E-12, RHOCNTRL = 1.7E-12, BETA = -45 
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6.0 Conclusions and Recommendations 


Two viable new techniques of momentum management for the the COG Planar 
Space Station configuration have been developed. These are the use of the 
core body for gravity gradient desaturation and the Implementation of four 
flat plate appendages at the end of long booms which would be used for aer- 
odynamic desaturation. The use of the core body, whose attitude is nomi- 
nally local vertical, for gravity gradient desaturation about the POP di- 
rection is a technique previously considered. The control law developed in 
this study shows that the core body can also be used to autonomously center 
POP momentum therefore minimizing momentum storage requirements. 

Using the core body to desaturate momentum TOP required a phased rotation 
about the vehicle X-axis to generated torque IOP with the desired magnitude 
and direction. All calculations using the core body for desaturation were 
autonomous and required only knowledge of the core body inertia matrix. 
With the exception of momentum sampling, calculations were only performed 
at orbital frequency and should Impose very little burden on a flight com- 
puter. Examination of the results showed that after equilibrium was at- 
tained, only small core body rotations at orbital frequency were required 
to control the momentum. 

As even these small rotations could interfere with payload pointing, an al- 
ternative method of momentum management using the core was developed. This 
method allowed the absolute momentum stored by the CMGs for an approximate- 
ly trimmed vehicle to slowly drift to a value inside the maximum storage 
sphere of four Skylab type CMGs. At that time, the control law would ini- 
tiate a maneuver, using the remaining CMG momentum, to remove the POP bias 
and re-center the POP momentum. The law would drive the IOP momentum along 
a diagonal to the opposite value of the IOP momentum circle. For these 
larger maneuvers, the profile was adjusted for finite vehicle maneuvering 
capability. This concept worked extremely well as maneuvers were required 
only about every thirtieth orbit. However the simulation did not include 
time varying external parameters such as e-angles and density functions. 
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To allow a trimmed vehicle to retain an inertial attitude for even longer 
periods of time, i f w aerodynamic desaturation appendages considered could 
be used. They do represent additional hardware requirements but an exami- 
nation of the results shows that the method is competitive with the bench*, 
mark, magnetic desaturation. The total weight would probably be less than 
the 170kg of the Space Telescope magnetic desaturation system and the power 
required much less than the 120 watts required at maximum magnetic mo- 
ment. However the control law utilizing these appendages does require at 
least an approximate value of the average air density during an orbit. 
Doubt exists as to the availability of this parameter but the control law 
developed has been shown to be relatively insensitive to atmospheric densi- 
ty uncertainties. One possibility is to autonomously obtain an estimate of 

p avg from the pop axis tr ^ m an 9 le * 

During periods of high atmospheric density, appendage desaturation is capa- 
ble of much larger torques than four Space Telescope torquer bars. During 
periods of very low atmospheric density, a higher degree of attitude trim 
might be required for the POP axis. I0P desaturation could be difficult 
since the appendages do not generate sufficient torque to counteract solar 
array GG bias torques at 45 degree e-angle. Appendages do not present the 
problems of magnetic contamination of payloads but potential problems could 
exist with 1 i ne-of -sight of payloads and attitude control position sen- 
sors. Due to the light weight and stiff construction proposed (Kevlar 
Boom), no attitude control structural mode induced control problems are an- 
ticipated and the transmission of torques to the core body should be quite 
rigid. By means of transformations stored in the attitude control comput- 
er, magnetic desaturation performs Independently of vehicle attitude since 
a magnetic moment can be defined in orbit inertial space. The appendage 
control law defined in this study only operates at the CDG Planar Space 
Station baseline attitude. Should appendage desaturation be desired at a 
variety of attitudes, the modified concepts could still be feasible but the 
appendages would require at least two-degrees-of-freedom. Additional 
studies could define the appendage/boom/actuator systems in greater detail 
and also evaluate control law performance over changing 8-angles and atmos- 
pheric densities. 
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